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Recent progress in the understanding of the nucleon is presented. The unpolarised structure functions are obtained 
with unprecedented precision from the combined Hf and ZEUS data and are used to extract proton parton distribution 
functions via NLO QCD fits. The obtained parametrisation displays an improved precision, in particular at low 
Bjorkcn x, and leads to precise predictions of cross sections for LHC phenomena. Recent data from proton-antiproton 
collisions at Tevatron indicate further precise constraints at large Bjorken x. The flavour content of the proton is 
further studied using final states with charm and beauty in DIS ep and pp collisions. Data from polarised DIS 
or proton-proton collisions are used to test the spin structure of the proton and to constrain the polarised parton 
distributions. 



1. Introduction 

The progress in understanding the structure of baryonic matter has been driven by the continuous improvement 
in the experimental ability to probe the matter with an increasing precision. Particle collisions in which one of the 
particles is preserved after the interaction is one of the most precise methods. For each collision the quantity of 
momentum transfered from the probe to the composed particle can be measured. This momentum, denoted by Q, 
reflects the virtuality of the exchanged boson and corresponds to a wavelegth given by A[fm] = 0.2/Q[MeV]. It 
gives the order of magnitude of the resolution with which the composed particle (e.g. the nucleon) is scrutinised 
in the collision - the higher the Q the finer the details probed in the interaction. A representation of the historical 
developments associated with the search for matter structure, from the discovery of the atomic nucleus in the classical 
Geiger-Marsden- Rutherford experiment, to the unique electron-proton collider HERA, is shown in figure 1 (left). For 
a momentum transfer of around 1 GeV the probed distances are comparable with the size of the proton. The inner 
structure of the proton was revealed in the first deep-inelastic scattering experiments at SLAC, 40 years ago. The 
"partons" were found and later identified with quarks and gluons within the framework of quantum chromodynamics 
(QCD). The "state of the art" is defined by the unique electron-proton collider HERA, where the centre-of-mass 
energy of 319 GeV allows the investigation of the proton with resolutions corresponding to 10" 18 m. 

The status of the proton structure investigations can also be represented in the plane spaned by Q 2 and the variable 
x, representing the fraction of the proton momentum carried by the interacting parton, as in figure 1 (right). The 
unique reach of HERA allows for the first time the investigation in the perturbative domain of the proton structure 
at very low x, down to 10~ 5 . This regime is particularly important for LHC physics, since the production of particles 
with moderate masses (a few hundred GeV's) may involve partons at low x. The measurements are complemented 
by the proton-antiproton collider Tevatron, where the investigations of W boson production constrain the proton 
structure at high x. Finally, the picture of the structure of the proton is also revealed in spin studies, performed using 
collisions of polarised beams. These studies aim to measure the various components contributing to the proton's spin 
and also provide independent tests of the strong interaction mechanisms. 



2. Proton structure and partons distribution functions 

The proton can be described in the infinite momentum frame as a collection of charged constituents (partons) that 
are "frozen" during the interaction time. In other words, in lepton-proton collisions, the time scale of the interaction 
between the incoming lepton and the proton's charged constituents is much shorter than the typical re-arrangement 
time of partons within the proton. Each interaction is characterised by the negative momentum transferred squared 
Q 2 and by the Bjorken x variable, defined as the fraction of the proton momentum carried by the interacting parton. 
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Figure 1: Left: The historical developments of the matter investigations. Right: The kinematic plane x-Q 2 , with the regions 
covered by the fixed target experiments, and the HERA and Tevatron colliders. The domain probed at LHC is also indicated. 

The understanding of high energy collisions involving protons within the QCD-improved quark parton model 
relies on the factorisation theorem, which states that the cross section is a convolution of a perturbative part, which 
describes the lepton-parton interaction at high energy, with a nonperturbative part that describes the momentum 
distributions of the proton inside the proton [1]. In the DIS case: 

<Jdis(x,Q 2 ) =C DIS ® f{x,Q 2 ) (1) 

The coeficient functions Cdis & re calculable at any order in QCD, while the parton distribution functions f(x,Q 2 ) 
(PDFs) encapsulate the parton content of the proton and are non-perturbative by definition. Their form cannot be 
calculated within QCD, nevertheless the evolution in Q 2 (i.e. their derivatives as a function of Q 2 ) are predicted 
by QCD via the DGLAP equations [2]. If the PDFs are parametrised as a function of a; at a given Q 2 = Qq, the 
evolution equations allow a recalculation at any Q 2 and therefore a prediction for any (x, Q 2 ) point of the DIS cross 
section that depends only on the PDF x-parameterisation at Qq. 

The factorisation ansatz is usually also applied, though not completely demonstrated, for proton-proton or proton- 
antiproton collisions. In this case, the perturbative coefficient functions, describing the partonic interactions (between 
quarks and gluons) at high energies, are convoluted twice with the PDFs, corresponding to the two interacting 
(anti)protons. The interaction scale (the equivalent of Q 2 in DIS) can be given by the measurement of the final state. 
For instance, the cross section for W boson production in pp collisions at Tevatron is expressed as: 

Vpp^WX = J dx 1 dx 2 C p artons->W <S> f{xi,M^) <g> f(x 2 , M$y) (2) 

The level of experimental precision achieved, together with a precise calculation of the perturbative coefficients allow 
these measurements to further constrain the PDFs in the domain at large x, as will be described below. 

Measurements in a large domain in x and Q 2 , as indicated in figure 1 (right), from DIS ep in fixed target and 
collisions at HERA, and from pp collisions at Tevatron are used to constrain the PDF parametrisation in a global 
fit. Each measurement at a given point (x,Q 2 ) is compared with the prediction obtained from a parameterisation of 
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the PDFs at a given scale Qq, then evolved at the actual Q 2 using the DGLAP equations, and combined with the 
perturbative coefficient functions. 



3. DIS cross sections and PDFs 

The kinematics of the lepton-nucleon reaction is conveniently described by the variables x, Q 2 , already introduced, 
and y, the inelasticity, representing the energy transfer between the lepton and the hadron system. The inelasticity 
y is directly related to the lepton-parton scattering angle, and can be used to study the helicity structure of the 
interaction. The kinematic variables are related by Q 2 = sxy, where s is the lepton-proton centre-of-mass energy. 

Neutral current (NC) and charged current (CC) processes are measured most recently by the HI and ZEUS 
experiments in e ± p collisions at HERA. The NC events contain a prominent electron and a jet of particles measured 
in the calorimeter, while in CC events only the jet is visible since the outgoing neutrino in not detected. 

Since a large domain in x and Q 2 is accessed, the NC cross section become sensitive to weak effects, besides the 
photon exchange, which dominates the cross section. The Z° boson exchange can be incorporated into the so-called 
generalised structure functions. The cross section is parameterised as: 

The helicity dependence of the electroweak interactions is given by the terms Y± = 1 ± (1 — y 2 ). The generalised 
structure functions F 2 and xF 3 can be further decomposed [3] as: 



K Q T?1 Z i f„.2 i „2\ f K Q \ p z 



with k 1 = 4-^^(1 — -rjjf ) in the on-mass-shell scheme. The quantities v e and a e arc the vector and axial-vector 

M z 1VI z 

weak couplings of the electron or positron to the Z°. The electromagnetic structure function F 2 originates from 
photon exchange only and dominates in most of the accessible phase space. 

The functions F 2 and xF 3 are the contributions to F 2 and xF% from Z° exchange and the functions F^ Z and xF^ z 
are the contributions from 7Z interference. These contributions are significant only at high Q 2 . For longitudinally 
unpolarised lepton beams the F 2 contribution is the same for e~ and for e+ scattering, while the xF 3 contribution 
changes sign as can be seen in eq. 3. The longitudinal structure function may be decomposed in a manner similar 
to F 2 . Its contribution is significant only at high y. 

In the quark parton model (QPM) the structure functions F 2 , F^ 2 and F 2 Z are related to the sum of the quark 
and anti-quark momentum distributions, xq(x,Q 2 ) and xq(x,Q 2 ): 

[F 2 ,Ff, F 2 \ = xYy v 2e q v q ,v 2 q + a 2 q ]{q + q} (4) 

and the structure functions xF^ z and xF% to the difference, which determines the valence quark distributions, 
xq v (x, Q 2 ): 

[xF^ z ,xF^] = 2xY^[e q a q ,v q a q ]{q- q} = 2x [e q a g , v q a q ]q v . (5) 

q q—u,d 

In equations 4 and 5, v q and a q are the vector and axial- vector weak coupling constants of the quarks to the Z°, 
respectively. 
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Figure 2: Left: the neutral current and charged current cross section as a function of Q 2 measured in e p collisions at HERA. 
Right: the charged current reduced cross section acc as a function of x for various Q 2 values measured in electron-proton 
collisions with different electron beam polarisations. 



The charged current (CC) interactions, e ± p -t'^X, are mediated by the exchange of a W boson in the t channel. 
The cross section is parametrised as: 

d 2 a cc (e ± p) G% 
dxdQ 2 2ttx 

with a± c (x,Q 2 ) = ±[Y + W±(x,Q 2 )TY_xW±(x,Q 2 )-y 2 w£ (x,Q 2 )] . 

The reduced cross section is denoted by a is the reduced cross section, Gf is the Fermi constant, Mw> the mass of 
the W boson, and W 2 , xW 3 and Wl, CC structure functions. In the QPM the structure functions and xW^ 1 
may be interpreted as lepton charge dependent sums and differences of quark and anti-quark distributions: 

W} = x(U + D), xW^=x(D-U), Wz=x{U + D), xW% = x(U — T>) , (6) 

whereas = 0. The terms xU, xD, xU and xD are defined as the sum of up-type, of down-type and of their anti- 
quark-type distributions, i.e. below the b quark mass threshold: xU — x(u + c), xD = x(d + s), xU — x(u + c), 
and xD = x(d + s) . 

The differential NC and CC cross sections as a function of Q 2 are shown in figure 2 (left) for e~p and e + p collisions 
at HERA. At low Q 2 the NC cross section, driven by the electromagnetic interaction, is two orders of magnitude 
larger than the CC cross section and corresponds to a pure weak interaction. At large Q 2 ~ z the two cross 
sections are similar. The largest Q 2 measurement corresponds to a resolution of S ~ 10~ 18 m, i.e. 1/1000 of the 
proton size. The agreement between the measurement and the prediction based on QCD improved parton model 
suggests no evidence for quark substructure. 

The double differential reduced cross section acc(x, Q 2 ) is shown in figure 2 (right). Two data sets are shown, 
with different polarisations of the electron beam. The enhancement of the cross section as a function of the beam 
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HERA I e + p Neutral Current Scattering - HI and ZEUS 

1 HI and ZEUS Combined PDF Fit 




Figure 3: Left: Examples of cross section combinations for three bins in i as a function of Q 2 . The precision of the combined 
experimental points is dramatically improved versus the individual measurements. Right: The combined HERA I data set 
compared with the HERAPDF 0.1 fit. 

polarisation <7p ±p = (1 ± P e )(J nopo i is visible. In addition, the CC processes are sensitive to individual quark flavours, 
which is especially visible at large Q 2 : the e + p collisions probe the d(x) quark distribution, while e~p are more 
sensitive to the u(x). This is a very useful feature of the CC processes compared to the NC, where the flavour 
separation is weaker and occurs only at high Q 2 where the Z boson exchange is important. 

4. H1 and ZEUS cross section combinations 

The cross section measurements are generally affected by statistical and systematic uncertainties. The latter 
can also be correlated among various measurements, leading to common shifts or coherent change of shapes in 
the determined PDFs. These uncertainties and their correlations should be properly taken into account in the fit 
procedures, in order the propagate the experimental uncertainties of the PDFs determinations into predictions for 
cross sections, for instance at LHC. The fit procedure effectively combines the information of the same measurements 
(for instance the cross section in the same phase space or bin). However, this combination is done via the regular 
assumptions of such fits: input parameterisations, tolerances among different data sets and so on. Alternatively, the 
measurements of the same observables can be combined in model independent way using a coherent procedure [4] that 
takes into account the statistical and sytematic errors. The main advantage of such procedure is the proper treatment 
of the correlated and uncorrelated errors. The obtained data set therefore reflects the experimental knowledge in the 
given measurement phase space and is factorised from the fit procedures, which can be therefore better understood 
and optimised. 

The average procedure is based on the minimisation of a x 2 function, whose free parameters are not only the average 
values in each measurement point but also the correlated systematics shifts from each data set. The minimisation 
procedure is designed to find the best estimation of the local average measurements, also taking into account also 
long range systematic shifts that may occur due to the correlated systematic errors of the measurements. This 
averaging procedure achieve a significant reduction of the long-range correlated systematic errors, as illustrated in 
figure 3 for three sets of measurements at low values of x. The gain in precision is significant and better than the 
simple combination of errors in each point. This can be explained by the fact that for each measurement point a 
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Figure 4: HERAPDF0.1 fit compared with MSTW and CTEQ fits. 

given correlated uncertainty of one set of measurements is constrained by the independent measurements from other 
data sets. This cross-calibration leads to a precise and coherent HERA data set to be used to determine the PDFs. 

5. QCD fits 

As explained above, a variety of data from DIS and hadronic collisions can be used to constrain the PDFs in the 
so-called global fits. The approach is based on the factorisation assumption and provides PDF parameterisations that 
can be used for further predictions, for instance at LHC. The consistency of this procedure provides also a stringent 
test of the QCD. Recent progress has been reported in including Tevatron data and refining the error evaluations [5] 
and in detailed studies of the precision of the observables and its correlation with the PDFs [6] . Such approaches are 
the main source for PDFs parameterisations since they incorporate all available information on the proton structure. 

However, it is well known that the global fits encouter and reveal tensions among various experimental data sets 
and, in absence of further experimental information, approximate statistical methods are used in order to cope with 
these inconsistencies. The most striking issue in the global fit procedures is the estimation of PDF errors, which 
has to be done with an somewhat arbitrary prescription of A\ 2 — 50 4- 100 in order to accomodate systematic 
differences between various data sets. The PDF parametrisations obtained in the global fits using various approaches 
are in good agreement most of the time, but there are also some cases where the predictions differ by more than the 
quoted PDF uncertainty [7, 8]. In addition, the various improvements of the treatment of the theoretical input have 
lead in the past to variations of the predicted cross sections at LHC by significantly more than the previously cited 
uncertainties [9]. This situation raises the question of the theoretical errors associated with the extracted PDFs. 

In these conditions, a useful approach is to use one coherent data set covering a significant area of the phase space, 
in order to study the QCD fit behavior, to compare to other data sets and to calculate predictions for observables 
at LHC. The newly obtained HERA combined data set has the advantage of very small systematic errors, which in 
fact no longer dominate. This makes the QCD fit procedure more robust and allows a test of the influence of the 
theoretical assumptions. The HERA combined preliminary data, described above, has been used by the HI and ZEUS 
structure functions working group to perform a stand alone NLO QCD fit, HERAPDF 0.1 [10]. The consistency of 
the input data set simplifies the determination of the experimental uncertainties on the PDFs, performed using the 
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Figure 5: The improvement in precision obtained using HERA data: prior to HERA data (left), using one HERA experiment 
(center) and using HERAPDF 0.1 parametrisation (right). Note that pre-HERA data do not constrain W production at LHC 
and the precision shown in the left figure fully relies on parameterisation extrapolation from large x measurements from fixed 
target experiments. 



more standard tolerance condition A% 2 = 1. 

In figure 3 (right), the NC DIS combined reduced cross section is compared to the HERAPDF 0.1 preliminary fit. 
The impressive precision of the combined data is reflected in the precision of the fit, which, extrapolated to lower 
Q 2 , is found to be in very good agreement with the fixed target DIS data. The precision of the new data lead to 
precise PDFs, in particular at low x, where most of the cross-calibration that occurs during the averaging procedure 
dramatically improves the systematic errors of the combined cross section. The obtained parameterisations together 
with the estimated errors are compared in figure 4 to the global fits based on separate data sets. 

Although the observed improvement in precision is mostly due to the improved combined data, the fit procedure 
itself, which is similar to the global fits, contains an intrinsic artificial constraint of the estimated error due to the 
functional form assumed for the PDFs at the initial scale Qq. Recent progress in the field of PDF parametrisation 
opens the way to an model independent assesement of the PDFs [11]. 

The parametrisation HERAPDF 0.1 is used to predict the W boson production in pp collisions at LHC. The 
improvement of the precision of this observable is shown in figure 5. The figure [12] displays the unique role of HERA 
for LHC physics and the significant gain in precision obtained recently from the HI and ZEUS combined data. More 
HERA data, in particular at high x and Q 2 and from high statistics jet production analysis, are expected in the near 
future. They will further improve the precision of PDFs and are therefore an essential input for the LHC physics. 



6. The direct measurement of the structure function F L at low x 

The longitudinal structure function Fl is a fundamental form factor of the proton, as can be seen from equation 3. 
Its contribution to the cross section is usually small, and only visible at large y where the weighting factor approaches 
the unity. Nevertheless, the structure function Fl is of particular importance. In the naive quark-parton model, the 
longitudinal structure function Fl = F 2 — 2xFi = (the Callan-Gross relation) reflecting the fact that longitudinal 
photons with helicity ±1 cannot absorb spin 1/2 quarks. This is possible only if a gluon radiation allows a spin flip 
of the interacting quarks. It can be shown that Fl is directly related to the gluon density in the proton [13, 14] 
xg(x) = 1.8[^-Fl(0Ax) — £2(0. 8x)] ~ —Fl meaning that at low x, to a good approximation Fl is a direct measure 
for the gluon distribution. 

A direct measurement of Fl can be performed if the cross section a(E p ) ~ F2(x, Q 2 ) + f(y) Fl(x, Q 2 ) is measured 
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Figure 6: The measurement of the longitudinal structure function Fl at HERA. 



at fixed x and Q 2 but variable y. This can only be performed if the collision energy ^fs is varied, for instance by 
reducing the proton beam energy at HERA. Then Fl(x,Q 2 ) can be directly measured with reduced uncertainties 
from the difference of cross sections: Fl ~ C(y) * (cr(Ep) — a(E 2 )). The data taking at lower proton energies took 
place at the end of the HERA run in 2007. The data collected at proton energies of 460, 575 and 920 GeV have been 
used to perform the first direct measurement of Fl in the low x regime. The first measurements in the medium Q 2 
range from 12 GeV 2 to 90 GeV 2 were published recently by the HI collaboration [15]. Preliminary measurements 
have been performed by HI [16] in an extended range of Q 2 up to 800 GeV 2 . The HI measurements averaged in x 
for each bin in Q 2 are shown in figure 6. The ZEUS collaboration presented measurements in the Q 2 range from 
24 to 110 GeV 2 [17]. The measured Fl is significantly different from at lowest Q 2 and are in agreement with 
the theoretical predictions based on various PDF parameterisations. It should be noted that in the presented phase 
space the gluon density and therefore the Fl predictions are well (but indirectly) constrained by other (inclusive) 
measurements. The measurements at lower Q 2 < 12 GeV 2 are experimentally challenging, but will be explored soon 
at HERA. They will shed light on various theoretical approaches predicting different values for Fl at low x [18] and 
together with other measurement in this kinematic regime will constraint even further the gluon component of the 
proton. 



7. Constraints from measurements in pp collisions at Tevatron 

Precise measurements of cross sections in pp collisions for processes well understood theoretically, can also be used 
to constrain the PDFs, in particular in the high x regime. Two such recent measurements are described here: the 
weak boson and jet production. 

Events with jets at Tevatron are produced by colliding quarks and gluons with high values of x. A recent mea- 
surement of the inclusive jet production was performed using a large statistics sample of CDF run II data [19], 
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CDF Run II Preliminary I L= 1 . 1 3 fb 1 




Figure 7: Left: Jet production cross section. Right: W production asymmetry. 



corresponding to 1.1 fb -1 , using the cone-based midpoint jet clustering algorithm in the rapidity region of \y\ < 2.1 
for jets with transverse momenta in the range from 60 GeV to 700 GeV. A similar measurement, using a kx jet clus- 
tering algorithm, was published by the DO collaboration [20] for jet rapidities in the range -2.4 to 2.4 in the transverse 
momentum range from 50 GeV to 600 GeV. The results are consistent with next-to-leading-order perturbative QCD 
predictions based on recent parton distribution functions. The dominant experimental uncertainty is related to the 
hadronic energy scale. Cross sections are measured as a function of the jet transverse momentum and in various 
regions of rapidity. The measurements are in good agreement with the NLO QCD prediction, as can be seen in 
figure 7, which shows results from DO. The main theoretical uncertainties originate from the PDFs and exceed the 
experimental precision. 

The measurement of W boson production at Tevatron also provides a probe of the proton structure. bosons 
are produced in pp collisions predominantly by the fusion of u{d) quarks from the proton and d(u) quarks from the 
anti-proton. As a simple example of the kinematics probed in this reaction at Tevatron and to compare with LHC, it 
can be noted that in order to produce a W boson at rest by the fusion of an quark from the proton with an d quark 
from the antiproton, the two quarks should have x — 0.2, while in the same configuration, the W will be produced 
at rest at LHC by partons with x = 0.005. 

The W + or W~ rapidity distributions reflect the quark momentum distribution inside the proton. As the u quark 
tends to carry a higher fraction of the proton's momentum than the d quark, the TV + (IV _ ) is boosted, on average, in 
the proton (anti-proton) direction. An asymmetry can be defined such that the charge of the W boson is correlated 
to the proton or anti-proton beam directions. The W charge asymmetry is defined as 

A(y w ) = (da(W +> /dy w - da{W~) / dy w ) / (da{W + ) / dy w + da(W~)/dy w ), (7) 

where yw is the rapidity of the W bosons. In the leading-order parton model, A(yw) is given approximately by 

A(y w ) = [u(xl)d(x2) - d{xl)u(x2)]/[u(xl)d(x2) + d(xl)u(x2)} (8) 

where = xoe ±yW and xq — Mw/i/s. Since the W four- vector cannot be unambigously reconstructed, a 
reweighting procedure is performed for the two reconstruction solutions obtained in each event by imposing the W 
mass constraint. The measurements were performed recently by both DO and CDF experiments [21, 22], each using 
each more than 1 fb -1 of run II data . The result obtained by DO is shown in figure 7 (right). The lepton charge 
asymmetry in Z° boson production leads also, in a similar way, to constraints on PDFs and has been measured 
recently with good precision at Tevatron [23] . 
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Figure 8: The charm structure function (left) and the reduced cross section for beauty production (right) measurements 
at HERA compared to theorectical prediction. The precision of the data, in particular that of the new measurements, is 
comparable or better with the difference among the various predictions and may lead to improvements in the proton structure 
function determination. 



The data from Tevatron measurements of the inclusive jets and W and Z charge asymetries were used in the most 
recent MSTW global fit and lead to significant improvements, in particular for the gluon in the large x domain [5]. 



8. Charm and beauty from the proton 

Of particular importance for the search of particles that couple to the mass at the LHC, especially for the Higgs 
studies, is the understanding of the heavy flavour content of the proton. In DIS, the interactions initiated by heavy 
quarks can be tagged by the presence of heavy flavour signals in the final state. Such interactions occur when gluons 
in the proton fluctuate to a pair of heavy quarks, one of which interacts with the photon (boson-gluon fusion), and 
therefore the study of heavy flavour interactions in DIS may also constrain the gluon content of the proton. The cross 
section for charmed tagged events can be written in a similar way as the inclusive DIS cross section, as expressed 
in equation 3. The main contribution to the charm (beauty) cross section originates from the structure function F£ 
(F^). At leading order, the F% and F\ structure functions are proportional to the momentum distribution of c and 
b quarks inside the proton. 

Measurements of heavy flavour production in DIS are performed at HERA [24] using two established experimental 
techniques: the tagging of D* mesons and the detection of particles emitted from long lived heavy hadrons. The 
measurements have been recently performed using the full available statistics from the HERA II run by the HI 
collaboration. The production of D* mesons has been mesured in the kinematic range of Q 2 from 5 GeV 2 to 
1000 GeV 2 [25]. The measurements are in good agreement with the calculations based on NLO QCD and with LO 
Monte Carlo programs. An extrapolation of the measured cross section from the visible phase space to the full phase 
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space is needed in order to extract the structure functions and to constrain the heavy quarks PDFs. The extrapolation 
is typically larger than a factor two and is dependent on theoretical assumptions. Therefore, a consistent theoretical 
framework has to be used in order to obtain the structure functions. The F 2 C structure function based on those 
precise measurements has been extracted recently [27]. 

The measurement of processes with heavy quarks in the final state based on their lifetime use precise silicon 
tracking devices in order to identify tracks with significant distance of closest approach to the beam axis, potentially 
originating from heavy quarks decays. A recent measurement by the HI collaboration [26] uses a multi-variate analysis 
and exploits the full data sets collected at HERA. A good agreement is obtained with the theoretical calculations, 
while the extrapolation factors to obtain the structure functions are closer to unity than in the case of D* method. 

A summary of the present status of the extraction of the heavy flavours structure functions F% and F| is presented 
in figure 8. The obtained precision is comparable to the theoretical calculation. The full potential of the HERA data 
will be obtained from the combination of HI and ZEUS final data and may play an important role for contraints on 
the gluon and on the heavy flavour content of the proton. 

Other measurements directly test the flavour aspects of the proton interactions. The measurement of events 
containing strange particles in the final state were performed at HERA [28, 29]. These measurements provide a 
new test of the QCD and allow a thorough investigations of strangeness production models. The link to the strange 
content of the proton is however difficult to establish since in the majority of events the strange particles are produced 
via the fragmentation process. Recent measurement at Tevatron [30] may also shed light on the strange content of 
the proton. The associated production of a W boson and a charm particle has been measured using soft muon 
tagging in a data sample corresponding to an integrated luminosity of 1.8 fb _1 . The measured process is described 
at leading order by the scattering of a gluon predominanlty with a strange quark. The measurement is found to be 
in good agreement with the theoretical prediction, but is so far statistically limited. A precision of about 15% on the 
strange quark content of the proton can be achieved with <~ 6 fb _1 of data, expected by the end of Tevatron run II. 



9. Study of the nucleon spin in polarised ep collisions 

The understanding of the nucleon structure cannot be complete without the understanding of the mechanism by 
which the quarks and the gluons compose the spin of the nucleon S z = 1/2. The various contributions to the nucleon 
spin can be expressed as following: 

S z = \ = + A.g(M 2 ) + Ll(^) + L 9 z (ji 2 ) . (9) 

Here AS (Ag) describes the integrated contribution of quark and anti-quark (gluon) helicities to the nucleon helicity 
and Lf (Lf) is the z component of the orbital angular momentum among all quarks (gluons) at a given scale /z 2 . 
The main puzzle has been the observation that, contrary to the naive expectation, the quark contribution does not 
account for the nucleon spin. 

The spin-dependent DIS cross section can be parametrised by two structure functions g\ and g2, where gi is 
negligible and g\ is given by: 

Here (e 2 ) — ^2 q e. 2 /N q is the average squared charge of all involved quark flavours, and Aq(x,Q 2 ) = q^(x,Q 2 ) — 
q^t(x, Q 2 ) is the quark helicity distribution for massless quarks of flavour q in a longitudinally polarised nucleon in 
the "infinite momentum frame" . 

The structure function gi is related directly to the cross section difference: cfll = ^( <7 ^ ~ <J ^)/'2 ; where longitu- 
dinally (i) polarised leptons (— >) scatter on longitudinally (X) polarised nuclear targets with polarisation direction 
either parallel or anti-parallel (=>, <=) to the spin direction of the beam. The relationship to spin structure functions 
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Figure 9: Left: The measured structure function gi for proton, deuteron and neutron. The precise measurement from HERMES 
is compared with measurements from other experiments. Center: The estimations of the gluon contributions Ag/g (from [32]). 
Right: Cross sections spin asymmetries for ty° and jet production in polarised pp collisions at RHIC, compared with a global 
fit. 
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(10) 



where 7 2 = Q 2 jv 2 , with v the energy of the virtual photon in the target rest frame. 

Measurements of <?i for the proton, deuteron and neutron are shown in figure 9(left) (from [31]). They can be used 
to extract the contribution of sea and valence quarks to the proton spin. Within some theoretical assumptions, this 
contribution is found to be AS(Q2 =5 GeV 2 ) — / dxgi(x,Q 2 ) ~ 0.33, which leaves a significant fraction for the gluon 
contribution to the proton spin. 

The gluon contribution to the proton spin can in principle be extracted, similarly to the unpolarised case, from the 
scaling violations of the polarised structure functions. Unfortunately, in the phase space accessed by the fixed target 
experiments - relatively large x and low Q 2 - scaling violations are expected to be small and the precision of the data 
limit this approach. An alternative method is based on the identification of hadrons in the final state with significant 
transverse momentum, enriching the selected sample in events initiated by direct gluon-photon interactions and 
thereby enhancing the sensitivity of the analysis to the gluon polarisation Ag/g. A summary of the various analyses 
leading to estimations of the gluon contributions Ag/g is presented in figure 9 (center, from [32]). The precision of 
these extractions are still limited. 

The spin of the proton is also investigated in polarised pp collisions at RHIC [33]. Recent measurements of the 
cross section spin asymmetries for processes with inclusive jets in the final state open the way to new constraints, 
in particular for the gluon polarisation. The correlation between the observed jets allows for an approximate recon- 
struction of the kinematics. The measurement of the spin asymmetry as a function of the jet transverse momentum 
allows the investigation of the integrated gluon contribution to the proton spin AG(Q 2 ) = J dxAg(x, Q 2 )- 

Global fits have been performed recently with the aim of determining the polarised PDFs [34] using measurements 
from polarised inclusive and semi- inclusive DIS and from polarised proton- proton collisions. The analysis shows 
that the RHIC data have a good sensitivity to Ag. The measurement of the double helicity asymmetry All, = 
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(c++ — er + _)/cr ++ — er H ) is compared with the fit in figure 9 (right). The fit indicates a positive Ag in a restricted 

region of x. The global fits are expected to also incorporate the analysis of final states with photons and W bosons. 
These analyses have a significant potential to enhance the sensitivity of polarised pp collision data to the gluon 
polarisation. In addition, the data collected from polarised pp collisions at higher energy ^/s^ = 500 GeV and lower 
energy ^/s^ = 64 GeV [35] will allow to explore a larger x domain, in particular at low x where Ag is unconstrained 
at present. The new measurements together with the expected refinements of the global fit procedures are expected 
to significantly improve the knowledge of polarised PDFs in the future. 

10. Outlook 

The research of nucleon structure is the scene of fast progress. Precise data from HERA constrain the parton 
distribution functions to an unprecedented level and provide solid support for predictions at the LHC in most of the 
phase space. The picture is complemented by Tevatron data at large x. The global fits exploit - and also prove - the 
universality of the parton distribution functions, combining the most precise results, the yield of which will increase 
in the years to come. Knowledge of the proton structure is profoundly enriched by spin studies in polarised collisions. 
Future proposed projects at higher energies like LHeC [38] and EIC [39], have the potential sharpen significantly the 
precision of nucleon structure investigations and to bring this fundamental research area into new paradigms. 
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